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Abstract. We have investigated the advection of a passive scalar quantity by 
incompressible helical turbulent flow in the frame of extended Kraichnan model. 
Turbulent fluctuations of velocity field are assumed to have the Gaussian statistics with 
zero mean and defined noise with finite time-correlation. Actual calculations have been 
done up to two-loop approximation in the frame of field-theoretic renormalization group 
approach. It turned out that space parity violation (helicity) of turbulent environment 
does not affect anomalous scaling which is peculiar attribute of corresponding model 
without helicity. However, stability of asymptotic regimes, where anomalous scaling 
takes place, strongly depends on the amount of helicity. Moreover, helicity gives rise 
^^ ■ to the turbulent diffusivity, which has been calculated in one-loop approximation. 
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1. Introduction 

During the last decade much attention has been paid to the inertial range of fully 
developed turbulence, which contains wave numbers larger then those that pump the 
energy into the system and smaller enough then those that are related to the dissipation 
processes pQ. Grounding of the inertial range turbulence have been created in the 
well known Kolmogorov-Obukhov (KO) phenomenological theory (see, e.g., PQI21E1)- 
One of the main problems in the modern theory of fully developed turbulence is 
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to verify the validity of the basic principles of KO theory and their consequences 
within the framework of a microscopic model. Recent experimental and theoretical 
studies indicate possible deviations from the celebrated Kolmogorov scaling exponents. 
The scaling behavior of the velocity fluctuations with exponents, which values are 
different from Kolmogorov ones, is called as anomalous and usually is associated with 
intermittency phenomenon. Roughly speaking, intermittency means that statistical 
properties (for example, correlation or structure functions of the turbulent velocity 
field) are dominated by rare spatiotemporal configurations, in which the regions with 
strong turbulent activity have exotic (fractal) geometry and are embedded into the vast 
regions with regular (laminar) flow. In the turbulence such phenomenon is believed 
to be related to the strong fluctuations of the energy flux which, therefore, leads to 
deviations from the predictions of the aforementioned KO theory. The deviations, 
referred to as "anomalous" or "non-dimensional" scaling, manifest themselves in singular 
(arguably power-like) dependence of correlation or structure functions on the distances 
and the integral (external) turbulence scale L. The corresponding exponents are 
certain nontrivial and nonlinear functions of the order of the correlation function, the 
phenomenon referred to as "multiscaling" . 

Although the theoretical description of the fluid turbulence on the basis of the "first 
principles", i.e., on the stochastic Navier-Stokes (NS) equation pQ remains essentially 
an open problem, considerable progress has been achieved in understanding simplified 
model systems that share some important properties with the real problem: shell models 
j3], stochastic Burgers equation |5. and passive advection by random "synthetic" velocity 
fields 0. 

The crucial role in these studies are played by the models of advected passive scalar 
field [7.. A simple model of a passive scalar quantity advected by a random Gaussian 
velocity field, white in time and self-similar in space (the latter property mimics some 
features of a real turbulent velocity ensemble), the so-called Kraichnan's rapid-change 
model [S], is an example. The interest to these models is based on two important facts: 
first, as were shown by both natural and numerical experimental investigations, the 
deviations from the predictions of the classical Kolmogorov-Obukhov phenomenological 
theory [U El Ej is even more strongly displayed for a passively advected scalar field than 
for the velocity field itself (see, e.g., 011111 HI] and references cited therein), and second, 
the problem of passive advection is much more easier to be consider from theoretical 
point of view. There, for the first time, the anomalous scaling was established on 
the basis of a microscopic model ^2]; and corresponding anomalous exponents was 
calculated within controlled approximations [T3| IT^j (see also review [Hj and references 
therein). 

The greatest stimulation to study the simple models of passive advection not only of 
scalar fields but also of vector fields (e.g., weak magnetic field) is related to the fact that 
even simplified models with given Gaussian statistics of so-called "synthetic" velocity 
field describes a lot of features of anomalous behavior of genuine turbulent transport 
of some quantities (as heat or mass) observed in experiments, see, e.g., [H] [IS] and 
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references cited therein. 

The term "anomalous scaling" reminds of the critical scaling in models of 
equilibrium phase transitions. In those, the field theoretic methods were successfully 
employed to establish the existence of self-similar (scaling) regimes and to construct 
regular perturbative calculational schemes (the famous e expansion and its relatives) 
for the corresponding exponents, scaling functions, ratios of amplitudes etc; see e.g. 
|16[ ITTj . Here and below, by "field theoretic methods" we mean diagrammatic and 
functional techniques, renormalization theory and renormalization group, composite 
operators, operator-product expansion and so on [17J. 

The feature specific to the theory of turbulence and simplified models associated 
with it is the existence in the corresponding field theoretical models of the composite 
operators with negative scaling (critical) dimensions. Such operators, termed 
"dangerous" in (1311312011211122], give rise to the anomalous scaling, i.e., the singular 
dependence on the infrared (IR) scale L with certain nonlinear anomalous exponents. 

Important advantages of the RG approach are its universality and calculational 
efficiency: a regular systematic perturbation expansion for the anomalous exponents was 
constructed, similar to the well-known e-expansion in the theory of phase transitions, 
and the exponents were calculated in the first order of expansion for passively advected 
vector fields (2H1 121] and in the second ^H] and third J2HJ orders of that expansion for 
scalar fields. Furthermore, the RG approach is not related only to the rapid-change 
model and can also be applied to the case with finite correlation time, anisotropy, the 
space parity violation, and, moreover, non-Gaussian advecting field [22]- 

The solution proceeds in two main stages. In the first stage, the multiplicative 
renormalizability of the corresponding field theoretic model is demonstrated and the 
differential RG equations for its correlation functions are obtained. The asymptotic 
behavior of the latter on their UV argument (r/l) (I is internal length) for r ^> £ 
and any fixed {r/L) (L is an outer length) is given by IR stable fixed points of those 
equations. It involves some "scaling functions" of the IR argument {r/L), whose 
form is not determined by the RG equations. In the second stage, their behavior at 
r C I is found from the operator product expansion within the framework of the 
general solution of the RG equations. There, the crucial role is played by the critical 
dimensions of various composite operators, which give rise to an infinite family of 
independent scaling exponents (and hence to multiscaling) . Of course, these both stages 
(and thus the phenomenon of multiscaling) have long been known in the RG theory of 
critical behavior. The distinguishing feature, specific to the models of turbulence is the 
existence of composite operators with aforementioned negative critical dimensions. Their 
contributions to the operator product expansion diverge at {r/L) — > 0. In the models 
of critical phenomena, nontrivial composite operators always have positive dimensions 
and determine only corrections (vanishing for {r/L) — > 0) to the leading terms (finite 
for {r/L) —> 0) in the scaling functions. 

Existence of regular perturbation schemes and accurate numerical simulations 
allows one to discuss, for the example of the rapid-change model and its descendants, 
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the issues that are interesting within the general context of fully developed turbulence: 
universality and saturation of anomalous exponents, effects of compressibility, 
anisotropy and pressure, persistence of the large-scale anisotropy and hierarchy of 
anisotropic contributions and so on. Moreover, it is interesting and important to study 
the helicity (violation of space parity) effects because many turbulence phenomena are 
directly influenced by them. 

In |22] the problem of a passive scalar advected by the Gaussian self-similar velocity 
field with finite correlation time J2S] was studied by the field theoretic RG method. 
There, the systematic study of the possible scaling regimes and anomalous behavior 
was present at one-loop level. The two-loop corrections to the anomalous exponents 
were obtained in [26J. It was shown that the anomalous exponents are nonuniversal 
as a result of their dependence on a dimensionless parameter, the ratio of the velocity 
correlation time, and turnover time of a scalar field. 

In what follows, we shall continue with the investigation of this model from the 
point of view of the influence of helicity on the scaling regimes and the anomalous 
exponents within two-loop approximation. 

Helicity is defined as the scalar product of velocity and vorticity and its non zero 
value expresses mirror symmetry breaking of turbulent flow. It plays significant role 
in the processes of magnetic field generation in electrically conductive fluid [23 I2H] 
and represents one of the most important characteristics of large-scale motions as well 
|291 I3U] . Despite of this fact the role of the helicity in hydrodynamical turbulence is not 
completely clarified up to now. 

The Navier-Stokes equations conserve kinetic energy and helicity in inviscid limit. 
Presence of two quadratic invariants leads to the possibility of appearance of double 
cascade. It means that cascades of energy and helicity take place in different ranges of 
wave numbers analogously to the two-dimensional turbulence and/or the helicity cascade 
appears concurrently to the energy one in the direction of small scales |HU 1321 ESJ- 
Particularly, helicity cascade is closely connected with the existence of exact relation 
between triple and double correlations of velocity known as "2/15" law analogously to 
the "4/5" Kolmogorov law [HUES]- Corresponding to [21] aforementioned scenarios of 
turbulent cascades differ each other by spectral scaling. Theoretical arguments given 
by Kraichnan [3E] and results of numerical calculations of Navier-Stokes equations [SZ] 
support the scenario of concurrent cascades. The appearance of helicity in turbulent 
system leads to constraint of non-linear cascade to the small scales. This phenomenon 
was firstly demonstrated by Kraichnan J3B] within the modelling problem of statistically 
equilibrium spectra and later in numerical experiments. 

Turbulent viscosity and diffusivity, which characterize influence of small-scale 
motions on heat and momentum transport, are basic quantities investigated in the 
theoretic and applied models. The constraint of direct energy cascade in helical 
turbulence has to be accompanied by decrease of turbulent viscosity. However, no 
influence of helicity on turbulent viscosity was found in some works [HHJ EHJ- Similar 
situation is observed for turbulent diffusivity in helical turbulence. Although the 
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modelling calculations demonstrate intensification of turbulent transfer in the presence 
of helicity 0U] direct calculation of diffusivity does not confirm this effect JJU 132] . 
Helicity is the pseudoscalar quantity hence it can be easily understood, that its influence 
appears only in quadratic and higher terms of perturbation theory or in the combination 
with another pseudoscalar quantities (e.g., large-scale helicity). Really, simultaneous 
consideration of memory effects and second order approximation indicate effective 
influence of helicity on turbulent viscosity [43J and turbulent diffusivity [HH EHl H3] 
already in the limit of small or infinite correlation time. 

Helicity, as we shall see below, does not affect known results in one-loop 
approximation and, therefore, it is necessary to turn to the second order (two-loop) 
approximation to be able to analyze possible consequences. It is also important to say 
that in the framework of classical Kraichnan model, i.e., the model of passive advection 
by the Gaussian velocity field with 5-like correlations in time, it is not possible to study 
the influence of the helicity because all potentially "helical" diagrams are identically 
equal to zero at all orders in the perturbation theory. In this sense, the investigation of 
the helicity in the present model can be consider as the first step to analyze the helicity 
in genuine turbulence. 

2. Field theoretic description of the model 

The advection of a passive scalar field 9{x) = 9(t, x) in helical turbulent environment is 
described by the stochastic equation 

d t e + v l d l = u o A6 + f, (1) 

where d t = d/dt, di = 8/8x4, ^0 is the molecular diffusivity coefficient (hereafter all 
parameters with a subscript denote bare parameters of unrenormalized theory; see 
below), A = d 2 is the Laplace operator, Vi = Vi(x) is the i-th. component of the 
divergence-free (owing to the incompressibility) velocity field v(x), and / = f(x) is 
an artificial the Gaussian random noise with zero mean and correlation function 

(f(x)f(x')} = 8(t - t')C{v/L), r = x - x', (2) 

where L denotes an integral (outer) scale. It maintains the steady-state of the system 
but the detailed form of the function C(r/L) is unessential in our consideration. In 
spite of the fact that in real problems the velocity field v(x) satisfies the Navier-Stokes 
equation, in what follows, we suppose that the statistics of velocity field is given in the 
form of Gaussian distribution with zero mean and correlator 

|^i^.(k) A^, k) exp[-^(t - if) + *k(x - x% (3) 

with 

D ^ U ' fe ) = T i 7 2-nW ' 1 M=^V ( 4 ) 

where k = |k|, D — g z/Q is a positive amplitude factor, g plays the role of the coupling 
constant of the model, an analog of the coupling constant X in the AqV? 4 model of 
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critical behavior fUJ Ej- in addition, g is a formal small parameter of the ordinary 
perturbation theory. The positive exponents e and rj (e = 0(r))) are small RG expansion 
parameters, the analogs of the parameter e = 4 — d in the Ao<^ 4 theory. Thus we have 
a kind of double expansion model in the e — rj plane around the origin e = rj = 0. The 
correlator (J3J is directly related to the energy spectrum via the frequency integral [22J 

2 

E(k) ~ A;^ 1 / dwD^fw, fc) ~ ^k 1 - 26 . (5) 

■/ MO 

Therefore, the coupling constant go and the exponent e describe the equal-time velocity 
correlator or, equivalently, energy spectrum. On the other hand, the constant u and 
the second exponent rj are related to the frequency uj ~ u UQk 2 ~ v which characterizes the 
mode k [46 . Thus, in our notation, the value e = 4/3 corresponds to the well-known 
Kolmogorov "five-thirds law" for the spatial statistics of velocity field, and rj = 4/3 
corresponds to the Kolmogorov frequency. For completeness, we remain ci-dependence 
in expressions © and (@J (d is the dimensionality of the x space), although, of course, 
when one investigates system with helicity the dimension of the x space must be strictly 
equal to three. To include helicity the transverse projector -P£(k) is taken in the form 

/§(k) = P tJ (k) + H Z] (k) = <% - kikj/k 2 + ipe^. (6) 

Here -Pij(k) = 8ij — kikj/k 2 represents non-helical part of the total transverse projector 
Pfj(\s). On the other hand, ifjj(k) = ipe^iki/k mimics the presence of helicity in the 
flow. Thus, formally, the transition to the helical fluid corresponds to the breaking of 
spatial parity, and, technically, this is expressed by the fact that the correlation function 
is specified in the form of mixture of a true tensor and a pseudotensor. The tensor e%n 
is Levi-Civita's completely antisymmetric tensor of rank 3 and the real parameter p, 
characterizes the amount of helicity. Due to the requirement of positive definiteness of 
the correlation function the absolute value of p must be in the interval \p\ £ (0, 1). Non- 
zero helical part proportional to p physically expresses existence of non-zero correlations 
(v • rot v). 

The general model (|3*|). (@J) contains two important special cases: rapid-change 
model limit when u — >• oo and g = go/^o = cons t, D v (u,k) — » g z/ A;~ d ~ 2e+r ', and 
quenched (time-independent or frozen) velocity field limit which is defined by uq — > 
and g ' = go/ito = const, D v (u,k) — > g , Qi , Q7c5(uj)k~ d+2 ~ 2£ , which is similar to the well- 
known models of the random walks in random environment with long range correlations; 
see, e.g., 03HSI- 

Using Martin-Siggia-Rose mechanism (IHj the stochastic problem (JTj)-(JH) can be 
treated as a field theory with action functional 

S(6, 9', v) = 9'D e 9'/2 + 9'[-d t + u A - (vA)}9 - vD; 1 v/2, (7) 

where 9' is an auxiliary scalar field, and Dg and D v are correlators (J2J) and (|HJ), 
respectively. In the action (JJJ) all the required integrations over x = (t, x) and 
summations over the vector indices are understood. The first four terms in Eq. (J2J) 
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Figure 1. (Left) Graphical representations of the needed propagators of the model. 
(Right) The triple (interaction) vertex of the model. Momentum k is entering into the 
vertex via field 9' . 

represent the Dominicis-Jansen type action for the stochastic problem (Q), (J2J) at fixed 
v, and the last term represents the Gaussian averaging over v. 

The model (J7J) corresponds to a standard Feynman diagrammatic technique with 
the bare propagators (99') q and (viVj)o (in the momentum- frequency representation) 

(%,k)gV,k)) = -iul^k* ' (vi(uMv j (u,k)) = Pf j (k)D v (u,k), (8) 

where D v (u>,k) is given directly by (jl|). In the Feynman diagrams, these propagators 
are represented by the lines which are shown in figure ^ (the end with a slash in the 
propagator (99') corresponds to the field 9', and the end without a slash corresponds to 
the field 9). The triple vertex (or interaction vertex) —9'vjdj9 = 9'vjVj9, where Vj = ikj 
(in the momentum-frequency representation), is presented in figure^ where momentum 
k is flowing into the vertex via the auxiliary field 9'. 

3. Renormalization group analysis 

The model (J7J) is logarithmic for e = rj = (the coupling constant go is dimensionless) 
and, in this case, possible ultraviolet (UV) divergences have the form of poles in various 
linear combinations of e and rj in the correlation functions. Using the standard analysis 
of quantum field theory one finds that all divergences can be removed by the only 
counterterm of the form 9'A9 22J. Thus, the model is multiplicatively renormalizable, 
which is expressed explicitly in the multiplicative renormalization of the parameters 
go, Mo, and u in the form 

i/ = vZ u , g = gp 2e+v Z g , m = up v Z u . (9) 

Here the dimensionless parameters g, u and v are the renormalized counterparts 
of the corresponding bare ones, /i is the renormalization mass (a scale setting 
parameter), an artefact of dimensional regularization. Newly introduced quantities 
Zi = Zi(g,u;d,p;e,T]) = Zi(g,u;d, p;s) are renormalization constants (note that if p 
is non-zero then d — 3) and, in general, contain poles of linear combinations of e 
and rj. However, as detailed analysis shows, to obtain all important quantities as the 
7-functions, /3-functions, coordinates of fixed points, and the critical dimensions, the 
knowledge of the renormalization constants for the special choice rj = is sufficient up 
to two-loop approximation (see details in [22]) ■ 
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The renormalized action functional has the following form 

S R (9, 9', v) = 9'D e 9'/2 + 9'[-d t + uZ x A - (vd)}9 - v J D" 1 v/2, (10) 

where the correlator D v is written in renormalized parameters. By comparison of the 
renormalized action ()10|) with definitions of the renormalization constants Zj, % = g, u, v 
(jHI one comes to the relations among them: 

Z v = Z±, Z g = Z u , Z u = Z v . (11) 

The second and third relations are consequences of the absence of the renormalization 
of the term with D v in renormalized action (|l()jl . The parameter p a the fields 9,9',v 
are not renormalized therefore Z p = Z$ = Zqi = Z v = 1. 

The issue of interest is, in particular, the behavior of the equal-time structure 
functions 

S n (r) = ([9(t, x) - 9(t, x')] n ), r = |r| = |x - x'| (12) 

in the inertial range, specified by the inequalities I << r << L (I is internal length). 
Here parentheses <> mean the functional average over the fields 9, 9', v with the weight 
exp(Sji). In the isotropic case, the odd functions S 2n +i vanish, while for S 2n simple 
dimensionality considerations give 

S2n(r) = v n r 2n R 2n (r/l,r/L,g 0} u ,p), (13) 

where R 2n are some functions of dimensionless variables. In principle, they can be 
calculated within the ordinary perturbation theory (i.e., as series in g ), but this is not 
useful for studying inertial-range behavior: the coefficients are singular in the limits 
r/l — > oo and/or r/L —>■ 0, which compensate the smallness of go, and in order to 
find correct infrared behavior we have to sum the entire series. The desired summation 
can be accomplished using the field theoretic renormalization group (RG) and operator 
product expansion (OPE); see (11112011221 for details. 

The RG analysis consists of two main stages. On the first stage, the multiplicative 
renormalizability of the model is demonstrated and the differential RG equations for its 
correlation (structure) functions are obtained. The asymptotic behavior of the functions 
like (fT2|) for r/l » 1 and any fixed r/L is given by IR stable fixed points g*,w* (see 
below) of the RG equations and has the form 

S 2 n(r) = u^ n r 2n (r/l)-^R 2n (r/L,p), r/l » 1 (14) 

with certain, as yet unknown, scaling functions R 2 n(r/L, p) = R 2n (l,r/L,g :¥ ,u*, p). In 
the theory of critical phenomena [TBI E! the quantity A[£2 n ] = — 2n + 7„ is termed the 
critical dimension, and the exponent j n , the difference between the critical dimension 
AfS^n] and the canonical dimension — 2n, is called the anomalous dimension. 

On the second stage, the small r/L behavior of the functions R 2n (r/L, p) is studied 
within the general representation (fT4*|) using the operator product expansion (OPE). It 
shows that, in the limit r/L — > 0, the functions R 2n (r/L,p) have the asymptotic forms 

R 2n (r/L)=Y,C F (r/L)(r/L)^, (15) 
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where Cp are coefficients regular in r/L. In general, the summation is implied over 
certain renormalized composite operators F with critical dimensions A n . In case under 
consideration the leading operators F have the form F n = (Vj#Vj#) n . 

We have performed the complete two-loop calculation of the critical dimensions of 
the composite operators F n for arbitrary values of n, d, u and p and obtain them in the 
following form: 

A = A (1) e + A ( V A (i) = -n(n-2)(d-l) 

" n + n ' " 2(d-l)(d+2) [ } 

where AW is critical dimension obtained in one-loop approximation. Interesting 
technical details of these two-loop calculations will be present elsewhere. 

Two-loop contribution A^ is rather cumbersome and can be found in |25j . 
The main and interesting result consists in the fact that although separated two- 
loop Feynman graphs of operators F n strongly depend on helicity parameter p, such 
dependence disappears in their sum, which gives rise to the critical dimensions A n . We 
can conclude that in two-loop approximation anomalous scaling with negative exponents 
(|16|) is not affected by the existence of non-zero helical correlations (vrotv) in turbulent 
incompressible flow. It turns out, however, that region of stability of possible asymptotic 
regimes governed by fixed points of RG equations, where anomalous scaling takes place, 
and effective diffusivity strongly depend on p. 

Let us analyze asymptotic regimes in detail. The structure functions and the other 
statistical averages of random fields 6, 9' satisfy linear differential RG equations with 
linear differential operator T>rg '. 

V RG = V fJL + /3 g (g, u)d g + (3 u {g, u)d u - 7„(g, u)V v . (17) 

Here D x = xd x stands for any variable x and the RG functions (the (3 and 7 functions) 
are given by well-known definitions and in our case, using the relations (jlljl for the 
renormalization constants, they acquire the following form 

lv =V^]nZ v , (18) 

f3 g =V IM g = g(-2e- V + 3 lu ), (19) 

f3 u = V^u = u{-n + lu ). (20) 

The renormalization constant Z v is determined by the requirement that response 
function G = (06') must be UV finite when is written in renormalized variables. In our 
case it means that it has no singularities in the limit e,rj — > 0. The response function 
G is related to the self-energy operator E, which is expressed via Feynman graphs, by 
the Dyson equation. In frequency-momentum representation it has the following form 

G(W, P) = — — ^^ r. (21) 

-iuj + uqP 1 - S(u;,p) 
Thus, Z v is found from the requirement that the UV divergences are canceled in (|2*Tjl 
after substitution uq = vZ v . This determines Z v up to an UV finite contribution, 
which is fixed by the choice of the renormalization scheme. In the MS scheme all the 
renormalization constants have the form: 1 + poles in e, f] and their linear combinations. 
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Figure 2. The one and two- loops contributions to the self-energy operator £. 

In contrast to the rapid-change model, where only one-loop diagram exists (it is related 
to the fact that all higher-order loop diagrams contain at least one closed loop which 
is built on by only retarded propagators, thus are automatically equal to zero), in the 
case with finite correlations in time of the velocity field, higher-order corrections are 
non-zero. In two- loop approximation the self-energy operator £ is defined by diagrams 
which are shown in figure El 

As was already mentioned, in our calculations we can put r\ = 0. This possibility 
essentially simplifies the evaluations of all quantities [221 EH] • 

Two-loop calculations of divergent parts of diagrams in figure 121 give the 
renormalization constant Z u and anomalous dimension / ~f u (|18jl in the form: 

Z v = 1 + 9 -A + 9 -B + 9 -C , lv = -2(0.4 + 2g 2 B). (22) 

see 1 

Here A = (1 — d)/(2d(l + u)) is the one- loop contribution to the constant Z v and 
anomalous dimension 7^ and the two-loop ones are 

(d-l)(d + u) 



B 



4d 2 (d + 2)(l + u] 



5 



C = - 



36(1 +u) 3 
(d-1) 2 1 
J 2 8(1 + u) 

ab„ 1 o(o+l)6(6+l) Ji 




^ 2^1--:-:-——!. (23) 



where 2-^1(0, b, c, z) = l+f^^+ "(c+iVi^ z +■ ■ ■ represents the hypergeometric function. 
We substitute d = 3 in the helical part (proportional to the p), but for completeness we 
remain the ^-dependence in the non-helical one. In addition we have introduced new 
notation g = gSd/(2u(2n) d ) (Sd = 2ix d l 2 /T(d/2) denotes the <i-dimensional sphere). 

From the expressions (|18|) - (|20|) and ()22j) we are able to find and clasify all fixed 
points #*,«* which satisfy equations: 

/3 9 {g*, u*) = p u (g*, «*) = 0. (24) 

To investigate the infrared stability of a fixed point it is enough to analyze the eigenvalues 
of the 2x2 matrix Q of first derivatives: Qij = df3 9i /dgj (gi = g,u). The anomalous 
scaling is governed by the infrared stable fixed points, i.e., those for which both 
eigenvalues are positive. 
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Classification and detailed analysis of all fixed points, determination of region of 
their stability and influence of helicity will be present elsewhere. Here we confine 
ourselves to the most interesting IR stable fixed point, where both parameters g*,u* 
acquire non-trivial values at rj = e: 

9, = ((9i 1) + (gi 2) +9i 3) P 2 )e)e, ff ?> = §(l + «.), 
8?= f H ! ? •2F 1 (l,l;l; 7r -^ ? ), (25) 




37T „ (1 1 5 , , 

«•"=-- •' F 'l.2-2 ! 2 ! (i+^l- (26) 

Actually, the equation (|23jl represents a line of fixed points in g — u plane. The 
competition between helical and non-helical terms appears which yields a nontrivial 
restriction for the fixed point values of variable u to have positive fixed values for variable 
9- 

4. Effective diffusivity 

One of the interesting object from the theoretical as well as experimental point of 
view is so-called effective diffusivity v. In this section let us briefly investigate the 
effective diffusivity u, which replaces initial molecular diffusivity vq in Q due to the 
interaction of the scalar field 9 with random velocity field v. Molecular diffusivity 
vq governs exponential dumping in time all fluctuations in the system in the lowest 
approximation, which is given by the propagator (response function) (jBJ). Analogously, 
the effective diffusivity v governs exponential dumping of all fluctuations described by 
full response function, which is defined by Dyson equation ()21j) . Its explicit expression 
can be obtained by the RG approach. In accordance with general rules of the RG 
(see, e.g., |H|) all principal parameters of the model go, «o and vq are replaced by their 
effective (running) counterparts, which satisfy RG flow equations 

dg n i- -\ du n , _. dv _ , , ,_,_, 

s~T = Pg{g,u), s— = f3 u {g,u) s— = -u^ u (g,u) (27) 

as as as 

with initial conditions g\ s =± = g, u\ s =i = u, P| a= i = v. Here s = k/fi, (3 and 7 functions 

are defined in (|T5|) -()20| ) and the running parameters g,u, and v clearly depend on 

variable s. Due to special form of (3- functions Q19j). (|2*0|) we are able to solve the last 

equation (|2*7)l analytically. Using the first equation (|27j) and (|TU|l one immediately 

rewrites the equation for effective diffusivity in the form 

dv 7 V dq , . 

— = 28 

v 2e + i] - 37^ g 

which can be easily integrated. Using initial conditions the solution acquires the form: 

-(^) 1/3 -(^) 1/3 - (») 

We emphasize that above solution is exact, i.e., the exponent 2e+r) is exact too. However, 
in infrared region k <^C // ~ l~ l , g — > g*, which can be calculated only pertubatively In 
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Figure 3. (Left) The dependence of i/* on the helicity parameter p for definite IR 
fixed point values u* of the parameter u. (Right) The dependence of v* on the IR fixed 
point u* for the concrete values of the helicity parameter p. The value p c = 4/V3- It 
is a special value related to the analysis of the stability of the scaling regime which is 
not discussed here. 



the two-loop approximation g* = 
gj in IJ29J1 we obtain: 

1/3 

v ~ v, I — — I k 



^g* £ / 



gPe 



■2c±ri 



P) 



+ (g* ' ■ !h 



V* = 1 



(3) 



£ 2 and after Taylor expansion of 



ii 2) + gi 3) )s 



3g 



(i) 



(30) 



Remind that for Kolmogorov values e — r\ — 4/3 the exponent in (|3U|) becomes equal to 
—4/3. Let us estimate the contribution of helicity to the effective diffusivity in the fixed 
point (|25j) . In this point e = r\ ((2s + 77) /3 = e) and 

(3 + m,) „ (, . 7 1 



V,; 



1-e 



np 



'30(1 + u 




12(1 + w*) ^ _± V 2 ' 2 ' 2 ' (l + «*' 9 ' (ol) 

In figure El the dependence of the v* on the helicity parameter p and the IR fixed 
point u* for Kolmogorov value of parameter e is shown . As one can see from these 
figures when w* — > 00 (the rapid change model limit) the two-loop corrections to v* = 1 
are vanishing. Such behavior is related to the fact that within the rapid change model 
there are no two and higher loop corrections at all. On the other hand, the largest 
two- loop corrections to the v* are given in the frozen velocity field limit (y* — > 0). 

Finally, let us analyze time-behaviour of retarded response function G =< 88' > in 
the limit t — > 00. 
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In frequency- wave vector representation G(u, p) satisfies Dyson equation (|2T|). Self- 
energy operator E is expressed via multi-loop Feynman graps and can be calculated 
perturbatively. We have found its divergent part up to two-loop approximation and 
calculated its finite part with the one-loop precision. 

Using Dyson equation we find response function in time- wave vector representation: 

G{t, P ) = / %-^G^ p) = / p , - e 7* . (32) 

J Z7T J Z7T —IUJ + 1/qP — L(Co>, p) 

In the lowest approximation E(u>, p) = 0, thus the integral can be easily calculated: 
Go(t, p) = 9(t)exp~ llJJrt . Here 6(t) denotes usual step function and u r is residuum in 
complex plain u in point —iv§p 2 . According to [S3], where analogical problem have been 
analyzed for turbulent viscosity, we suppose that this situation remains the same for 
the full response function G, i.e., the leading contribution to its asymptotic behavior for 
t — > oo is determined by the residuum u = u r , which corresponds to the smallest root 
of dispersion relation 

G _1 (u;, p) = — iuj r + uop 2 — £(toV, p) = 0. (33) 

It is advantageous to rewrite the last relation in dimensionless form: 

iu r . . E(u>, p) , , 

l-z-I(l,z) = 0, z = — ^, J(M,«7=-^fS 34 

is p 2 is p 2 

which after renormalization can be rewritten in the fixed point g* (125(1 as follows 

__ iu r 

-L Z* -** — U , Z* = — — , v^^J 

up 2 

where v is effective diffusivity ((30|) and /* = I*(l,z*,g*) is renormalized (finite) part of 
dimensionless self-energy operator / at the fixed point g*. 
Hence decay law G (t, p) ~ exp~ uop * is changed into 

G(t, p) ~ exp"^ 4 = exp' z * pp2t t ->■ oo . (36) 

To find the residuum u r (or, equivalently, z*) it is necessary to calculalate quantity I*. 
In one-loop (linear in e) approximation it can be written in the form: 

-i _ j_i 

(37) 



with 



X 



~ 9 *L 


(1 -x 2 ) — dxl 






/•OO 

I dk 




k 




0{k-l)~ 

(1 +M*)/C_ 


1 — z*+ (1 + u*)k 2 - 


-2kx 



(38 



Generally, the root z* can be complex and in one-loop approximation it has form 

Z* = z{ + 1Z2, z\ = 1 + X±S, Z* 2 = X2£ ■ (39) 

With our guaranteed precission I* is linear in e, therefore on the first sight it seems 
that in the last integral it is enought to take z* — 1 (<?* ~ e), but, actually, for its correct 
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calculation we need to remain imaginary part X2S. Then the integral (|38|) can be easily 
calculated by means of Sokhotsky's formula: 



lim = ^fin5(y) + P( — 

£->o+ y ± ie y' 



(40) 



where P(-) denotes the principal value of the integral. Integration over angle x gives 
final result for dimensionless self-energy operator I*: 

_ g m f V^Tjd + 1) (7 + j>{l + J ) + 2 In |1 + u*\) . d -l 

"■~"[i + 0\ d 2^(^)r(l + i) ±m 2d 







(41) 



where 7 is Euler's constant and ip(z) is digamma function defined as ip(z) = T'(z)/T(z). 
Successful calculation of integral / allows one to determine the residue z* (13"HJ) . 
Comparison of real and complex part of both sides of ()35j) gives the following terms in 
real space d = 3 and in the fixed point g* = 3(1 + u*)/2 (see (J2SJ) 



Xl= 3 +21n ^~ 



7T 
X9 = ± — . 

2 



(42) 



Due to the existence of two complex conjugate values z* the response function 
G(t,p 2 ) can be written in the asymptotic limit t — > 00 in the following final form 

G(t,p 2 ) = e~^f p2 ~ H sm(v fP 2 - £ t), 



(43) 



where 



v, 



eff 



1-e 



8 , n1 1 + M* (3 + M*) 

— him ; — 

3 2 30(1 + w*) 3 



1 1 5 



Vf 




2'2'2'(1 + ^) 2 
1/3 




2 \3(1 + u*)ey 
It's clear that the exponential damping is accompanied by the oscilations. 



5. Conclusion 

We have studied the advection of scalar field by turbulent flow in the frame of extended 
Kraichnan model and investigated the influence of helicity on anomalous scaling, 
stability of asymptotic regimes and effective diffusivity. Such investigation is useful 
for understanding of efficiency of simplified models to study the real turbulent motions 
by means of modern theoretical methods including renormalization group approach. 
Actually, we performed two-loop calculations of divergent parts of Feynman graphs, 
which are necessary to achieve multiplicative renormalization of equivalent field theoretic 
model. We have shown that anomalous scaling, which is typical for the Kraichnan 
model and its numerous extensions [2E1 OH] , is not violated by inclusion of helicity to 
the incompressible fluid. On the other hand, stability of asymptotic regimes, values 
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of fixed RG points and turbulent diffusivity strongly depend on amount of helicity It 
can be easily see from (J3~T|) that helicity enlarges turbulent diffusivity and high order 
contributions lead to the appearance of oscillations in response function ()32j) . 
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